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Abstract. Dynamical chiral symmetry breaking and confinement are two crucial features of Quan- 
' turn Chromodynamics responsible for the nature of the hadron spectrum. These phenomena, pre- 

C\ , sumably coincidental, can account for 98 % of the mass of our visible universe. In this set of lectures, 

I shall present an introductory review of them in the light of the Schwinger-Dyson equations. 
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5?: MOTIVATION 

Q_i| For millennia, mankind has always tried to explain how the world works. Every piece 

of knowledge, every new observation, every brilliant idea have modeled a masterpiece 
of that spot we like to call our home -the Universe. A common denominator of ancient 
cosmogonies has been to provide the elements necessary to explain what are the compo- 
nents of the Universe and how these ingredients sustain the delicate dance of the cosmos. 
Our modern understanding about the composition of the universe reveals that most of 
our home is dark, either in the form of dark energy, about a 74% of it, or as dark matter, 
another 22% [1]. Only a small portion of the Universe is what we see, and it is what has 
£Sj ■ delighted us since ancient times. 

Science, and particularly, physics, usually take a reductionist path to unveil the deep- 
est secrets of the Universe. In order to understand galaxies, we need to know everything 
about stars and planets: what are they made of, how they evolve in time, how they col- 
lapse. This, in turn, requires us to learn about molecules, atoms, nuclei, nucleons, and, 
ultimately, quarks. Our hope is that if we manage to explain the fundamental building 
blocks, that is, if we manage to know what are their basic properties, we should be able 
to explain what is the Universe made of. And we also need to know how these funda- 
mental blocks are glued together as to keep the Universe flowing. These are the basic 
questions that elementary particle physics tries to answer. 

Nature has given mankind the ability to ask questions. Most of the progress we have 
achieved was done answering questions, not necessarily wise man questions; the simpler, 
the better. One of such has to do with one property of matter very familiar to us in 
our everyday experience: Is mass a fundamental property of matter? At first sight, that 
questions is rather odd: We are so used to relate mass to matter, that we cannot think of 
one without the other. And it is true! Our everyday experience is ruled by the classical 
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laws of physics, elegantly summarized by Sir Isaac Newton [2]. In particular, Newton's 
second law establishes that a body accelerates proportionally to the force acting on it 
divided by its mass, in such a way that the motion of a body is intimately related to its 
mass. Furthermore, from Lavoisier we know that in every process, the mass is conserved. 
Thus, the Newtonian vision is that the mass of a body is the quantity of matter arising 
from its density and bulk conjointly. This brings us to the conclusion that the mass is, 
indeed, a fundamental property of matter at the classical level. 

It took one more genius to ask a simple question that crumbled the Newtonian vision 
of mass: "Is the energy of a body a measure of its energy content?" Answering this 
question, Albert Einstein [3] introduced perhaps the best known law of physics, 



where Eq is the body's rest energy. To be fair, that is not precisely the most famous form 
of the equation. That would be Eq = mc 2 , which somehow prejudices to think of energy 
as a mass, and not the other way around. The realization of Einstein's conception of 
mass takes place everyday in particle accelerators, where from the collisions of particles 
at sufficiently high energies, much heavier particles are produced. 

The difference between the Newtonian and the Einstenian conception of mass at the 
fundamental level is attainable only at the subatomic level (see, for instance, Ref. [4]). 
For example, consider the energy of a hydrogen atom in its ground state and think of it 
as being just its "Einstenian" mass. On the other hand, consider the mass of the proton 
and of the electron and add them in the "Newtonian" way. The difference of these two 
masses corresponds to the ionization energy of the hydrogen, which is a difference of 
i^(10~ 8 )! This tells us that at the atomic level, still the mass is Newtonian. A similar 
reasoning applies to the nucleus. Our empirical knowledge shows that most of the mass 
of a nucleus comes from the sum of the mass of its nucleons. It is hardly a 1% difference, 
for example, the difference between the mass of a nucleus that undergoes an a-particle 
emission as compared to the mass of the products of the decay. Such a difference is due 
to the nuclear binding energy. The mass of a nucleon, however, is a completely different 
story. Light quarks acquire, through the Higgs mechanism, current masses of, say, 3-5 
MeV, such that the sum of the masses of three light quarks is around 10 MeV, 98% less 
than the mass of the proton! This means that the mass of a nucleon is "Einstenian". 

A warning note is necessary here. It is usually quoted that the Higgs mechanism [5] 
(based upon the ideas of Spontaneous Symmetry Breaking of Nambu [6], whom was 
recently Laureated with the Nobel Prize of Physics 2008), is the origin of the mass in 
the Universe. This is far from being true. Although the Higgs mechanism accommo- 
dates the mass of the gauge bosons of the electroweak Standard Model respecting its 
gauge symmetry, and it also gives, say, electroweak masses to quarks and leptons, it cer- 
tainly cannot explain around 98% of the mass of our visible Universe, which is mostly 
composed of protons and nucleons. Understanding the Higgs mechanism is necesary to 
explain the existence of atoms and how do they form stable structures, a question as old 
as human kind itself. One of the physics goals of LHC [7] is precisely to find the Higgs 
boson and determine its propeties. Yet, the origin of mass of the visible Universe has an 
entirely different explanation. 



The origin of most of the mass we can see is explained by the strong interactions 
of QCD [8]. The story goes as follows: Our inability to observe colored excitations in 
a detector tells us that valence quarks are confined inside hadrons. Furthermore, when 
they are very close to each other, that is to say, when they are very energetic, they behave 
as free particles. This is the well known property of asymptotic freedom of QCD [9]. 
But when they are less energetic, namely, at a distance scale of the order of a nucleon 
size, their respective color clouds start to overlap. In these clouds, because of the strong 
interactions, virtual quarks and antiquarks start condensing, forming the chiral quark 
condensate. This condensate provides the valence quarks of the nucleon of a dense, 
sticky medium in which they propagate as though their masses were some 300 MeV. 
Adding the mass of three of these guys explains the value of the mass of a nucleon 
and hence of the visible Universe! The dynamically generated mass actually measures 
the confining energy of the quarks inside a nucleon. The phenomenon just described is 
known as the dynamical breaking of chiral symmetry, and occurs to light quarks even if 
their current masses vanish [10]. The heavier the quark, it becomes less sensitive to the 
dynamical effects of color interactions at low energies. Therefore, these heavy quarks 
acquire their masses through the Higgs mechanism instead. 

In this contribution I shall present a particular Framework to explore the details of the 
phenomena of confinement and dynamical chiral symmetry breaking, the Schwinger- 
Dyson equations [11], sometimes also referred to as Dyson-Schwinger equations. For 
this purpose, I shall work out a toy model of QCD which, however, is interesting by 
own merits, QED3, that is, the standard quantum electrodynamics restricted to a plane. 
Students are encouraged to follow the Approach and develop the necessary techniques 
to unveil the secrets of the origin of mass. 

FRAMEWORK 

A natural framework to study the strong interactions of QCD is, indeed, a lattice simula- 
tion. Lattice gauge theory has provided an accurate description of the hadron spectrum, 
within a 10% error of the experimental data [12], which confirms the physical picture of 
the origin of mass described above [8]. A different, and in many senses, complementary 
platform to study non perturbative phenomena in Quantum Field Theories is provided 
by the Schwinger-Dyson equations (SDEs) [11]. These are the field equations of a given 
theory, and form an infinite tower of relations among the Green's functions of the theory. 
In order to address the issues of confinement and dynamical chiral symmetry breaking 
(DCSB), and for the sake of simplicity, instead of solving the SDEs for QCD, I shall 
restrict myself to a simpler and more familiar theory, QED, starting from the Lagrangian 

Se = xfr(if(d^-ieA^) -mo) W - \f^ v - -j(^A") 2 , (2) 

where all the quantities carry their usual meaning. 




FIGURE 1. Perturbative Expansion of the Fermion Propagator 




FIGURE 2. Fermion Self Energy £(p) 

In order to understand the physical content of the SDEs (see, for instance [13, 15]), 
consider the perturbative expansion of the fermion propagator Sp(p) in QED, Fig. (1). 
There appear three types of radiative corrections, to the fermion propagator itself, to 
the photon propagator and to the fermion-photon vertex. Defining the self-energy as in 
Fig. (2), all these corrections are incorporated in this object, in such a fashion that the 
perturbative expansion in terms of £(p) now reads 

Sf(p) = Sf\p) +sf\ P n P )sP(p) +s£\p)Z(p)S<g\p)I,(p)sP(p) + ... (3) 



and is shown in Fig. (3). Here, SjP (p) is the bare fermion propagator. Factorizing Sp (p), 
the remaining expression is a geometric series in T,(p)Sp(p), which adds to 

S {0) (d) 

Sf( P ) = Sf (P 0) . (4) 

Equivalently, leaving the bare propagator apart and factorizing Sp (p)E(p)Sp (p), we 
again reach at a geometric series in Z(p)S ( p \p), and hence 

o(0)/ \ 

s F (p)=s?( P )+sP(pn P ) s ru> (5) 

l-Z(p)sP(p) 




FIGURE 3. Perturbative Expansion of the Fermion Propagator in terms of T.(p) 




FIGURE 4. Schwinger-Dyson equation for the Inverse Fermion Propagator in QED 



Comparing eqs. (4) and (5), we obtain 

S F (p) = S { °\p) +sP( P n P )S F (p) . (6) 

This is the SDE for the fermion propagator. On multiplying this expression by Sp 1 (p) 
from the left and Sp l (p) from the right, we arrive at the most familiar SDE for the inverse 
fermion propagator 

Sp 1 (p)=S^-\ P )-Z( P ), (7) 

and is depicted in Fig. (4). 



This corresponds to the expression 

Sp\p) =4 0) ~\p)~ie 2 J -^- 4 fS F (k)T v (k,p)A^(k-p) , (8) 

where e 2 is the electromagnetic coupling, A^y represents the complete photon propaga- 
tor and r v the full fermion-photon vertex. This is the equation we need to solve in order 
to study DCSB. 



Notice that the Green's functions involved in E(p) obey their own SDE. The SDEs for 
the photon propagator and the fermion-photon vertex are depicted in Figs. (5) and (6), 
respectively. Observe that the two-point Green's functions are coupled to one another 
and to a three-point function, which, in turn, is coupled to them and to a four-point 
function and so on. Hence, the set of SDEs, which are the field equations of QED, form 
an infinite tower of relations among the Green's functions of the theory. 

Now, we need to find out which kind of unknowns we are dealing with before we 
even try to find them. Let us notice that the complete fermion propagator, Sp(p), can 
be expressed in terms of two unknown scalar functions of the momentum squared, 



FIGURE 5. Schwinger-Dyson equation for the Photon Propagator in QED 





FIGURE 6. Schwinger-Dyson equation for the Fermion-Boson Vertex in QED 




multiplied each by the tensor structures 1 and p 1 which appear in the Dirac equation. 
For reasons that will become clear shortly, we write the propagator as 

and refer to F(p) as the fermion wavefunction renormalization and M(p) as the mass 
function. Notice that the pole of the propagator, i.e., the mass of the particle, is located 

at p 2 = M 2 (p). The bare propagator Sp(p) corresponds to F(p) = 1 and M(p) = niQ. 
The photon propagator, in turn, can be written as 

A**(«) = ^(«**-^)+^, GO) 

where G(q) represents the unknown photon wavefunction renormalization and % is the 
usual covariant gauge parameter. £ = labels the Landau gauge, and £ = 1 the Feynman 

gauge. The bare propagator A^l(q) is obtained setting G{q) = 1. Finally, the fermion- 
photon vertex can be expressed as 

r^(M = I>(Mv/\ (ii) 

!=1 

where Vi(k,p) are unknown scalar functions and V; are the basis vectors formed from 
the products of each the three vectors k* 1 , p^ and with each of the 4 spin structures 
1, p" and ftp". In summary, there are 15 unknown functions in eq. (8). Fortunately, the 
gauge principle tells us that not all of them are independent, as they are related through 
gauge identities like Ward-Green Takahashi identities [14] and so on. A recent review 
on different truncations of SDE can be found in [15]. 



Confinement and DCSB are usually studied from the fermion propagator. This means 
that if we want to truncate the tower of SDEs at the level of eq. (8), reasonable assump- 
tions have to be made on the form of the photon propagator and the fermion-photon 
vertex. Below we shall consider one of such truncation schemes which, on one hand, 
makes the analysis very neat and, on the other hand, shows the key features of QCD that 
explain the origin of mass. 

APPROACH 

In order to study DCSB and confinement, let us work in a toy model of QCD: QED3. 
That is, let us work with the ordinary QED, but restricted to a plane. Such a theory is 
super-renormalizable, so that we don't have to deal with regularization issues in a trun- 
cation of SDEs (see, for example, Ref. [15] for the implications of this fact). Contrary 
to the ordinary QED, the planar theory exhibits confinement [16], in a manner such that 
it resembles more to QCD that to QED itself. Even more, it has been argued very re- 
cently that both the restoration of chiral symmetry and the confinement/deconfinement 
phase transitions take place simultaneously in this theory [17]. Furthermore, to the cor- 
responding Lagrangian we can add a Chern-Simons term, which induces a topological, 
gauge invariant mass to the photons, making the structure of the theory even richer (a 
discussion of the most general QED3 Lagrangian can be found in Ref. [18]). QED3 has 
multiple applications not only because it emerges as the infinite temperature limit of 
QED [19], but there are many condense matter systems, ranging from high-T c . super- 
conductors [20] the quantum Hall effect [18, 21], and more recently, graphene [22], de- 
scribed by this theory. In the field of DCSB, QED3 has provided a popular battleground 
for lattice, SDEs and other studies, particularly regarding a possible critical number of 
fermion families for which the phenomenon ceases to take place [23]. 

We start our study of DCSB from the SDE for the quark propagator Sf(p), which in 
QED3 is given by 

SfHp) = 4 0) ~\p)~ie 2 J ^fS F (k)r v (k,p)A^(k-p) . (12) 

Here, e 2 is the dimensionful coupling of the theory that sets its natural scale. This means 
that any new mass scale, like the dynamical mass, should be proportional to e 2 . To 
emphasize the phenomenon of DCSB, we consider massless fermions to start with, 
setting mo = in the bare propagator. In order to truncate the tower of SDEs at this level, 

we need ansatze for T v and A^y Let us start by in neglecting fermion loops, A^ v — > A^° v . 
This is the quenched approximation of QED, which is pictorially represented in Fig. (7). 

With a suitable choice of the fermion-photon vertex, the SDE for the fermion propa- 
gator can be solved self-consistently. Possible the simplest interaction we can consider 
is the bare interaction, T v (k,p) = y v '. This is the so-called rainbow approximation, de- 
picted in Fig. (8), and under which the gap equation reads 



FIGURE 7. Quenched truncation of the Schwinger-Dyson equation for the Fermion Propagator in QED 
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FIGURE 8. Rainbow truncation of the Schwinger-Dyson equation for the Fermion Propagator in QED 
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Here we have set a = e 2 /(4n), as usual. This is a matrix equation which can be 
converted into a system of scalar integral equations for F(p) and M(p) after multiplying 
it by p and 1. However, in the Landau gauge F(p) = 1 (see, for example, Ref. [24]), 
so, in this gauge we are left with a single equation for M(p) alone. Such an equation is 
derived taking the trace to eq. (13) with F(p) = 1, yielding 

After contracting with the photon propagator, we carry out a Wick rotation to Euclidean 
space with the following prescriptions, 

(k°,p°) -> (ik° E ,ip° E ), 

{k\ p 2 , (k- P ) 2 ) - (-4, -p 2 , -(k- P ) 2 ) , 

d 3 k -> i [ d 3 k E . (15) 



For ease of notation, we shall drop the subscript E from every Euclidean space expres- 
sion. Hence, we are finally left with the following expression for the Euclidean gap 
equation in QED3, 

at d 3 k M(k) 

M{P) = T 2 ! (k- P ) 2 k 2 + M 2 (k) - (16) 

We can still go further making reasonable simplifying assumptions [25]. Since the 
DCSB is an infrared phenomenon, we expect the dynamically generated mass to serve 
as an IR cut-off for the gap equation. We can then linearize the above expression setting 
in the denominator of the integrand 



M 2 (k)=M 2 (0)=m 2 dyn , 



(17) 



which yields the linearized gap equation 

XMl N a f d 3 k M(k) 

Notice that M(p) = is a (trivial) solution of this equation, and would correspond to 
that derived in perturbation theory. However, we are interested in a non trivial solution, 
which we shall look for making use of both analytical and numerical techniques. 



Analytical treatment 

In order to solve eq. (18) analytically, define the function x(p) through [25] 

M( P ) = ( P 2 + mj yn )x(p). (19) 



dyn 

Then we are left with 

d 3 k 

Now, defining the Fourier transform of x(p)> 

d 3 p 

we observe that it satisfies the differential equation 



(p 2 + ™%n)x(p) = ^ 2 Jj j ^X(k). (20) 
sform of x(p)> 

cp(r) = J -^X(p)e ip \ (21) 



whose solution is 



<p(r) = Ce- md y" r . (23) 

On Fourier transforming back the result to momentum space, the normalization constant 
C is found demanding the infrared behavior of M(p — > 0) = m dyn . The resulting mass 



function is then 



3 

m. 



M(P) = . (24) 

and is depicted in Fig. (9). We observe that the mass function falls-off as M(p — > oo) ~ 
1 / p 2 . This is consistent with the expression for the chiral condensate, 

{f¥) = l im , (25) 

derived from the Operator Product Expansion (see, for example, [24]), since it demands 
the ultraviolet behavior M(p) ~ l/p 2 in order to render the condensate momentum 
independent. These limits will become relevant to test the validity of the numerical 
procedure that we are about to develop below. 
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FIGURE 9. Analytical solution to the linearized Schwinger-Dyson equation for the Fermion Propagator 
in QED3. The scale of the graph is set by m c [ yn = 1 



Numerical treatment 

Consider again eq. (18) and select a reference frame where the angle between mo- 
menta k and p is 6. Then we can write d 3 k = k 2 dksm0d6d(j) for G (0,2tt), 6 G (0, n) 
and k G (0, °°) . The integration over gives a factor of 2tz, and the integration over is 



sin0 



1 



/o k 2 + p 2 — 2kpcos6 2kp 
Gathering things together, we have 

kM(k) 



In 



k + p 



k — p 



„, v a f kM(k) 
M(p) = — / dk- — 4- In 
npJo k 2 + m 2 



k + p 



k — p 



(26) 



(27) 



Observe that there is an integrable singularity in the logarithm for k = p. In order to 
avoid numerical complications coming from this singularity, we can approximate 
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k + p 



k — p 



?n 2k 
±f#(k-p) + —#(p-k) 
k p 



(28) 



where &(x) is the step function. Such an approximation is valid for k~^> p and k <^ p, 
and yields the simplified SDE given by 

r dk ^k m - pH k - Hl ,- k) \ . ( 29) 

np Jo k 2 + mf lyn [k p J 

This is the equation we need to solve numerically. 

Set e 2 = 1 to start with. Then, select an integration interval (?c, X) such that 

M(p) * ^ f'dk-^-U »(k-p) + -*(p-k)\ . (30) 
JCpJk k 2 + m 2 yn [k p J 

This integral equation is of the form 

M(p)= dkKer(k,M(k);p) , (31) 

Jk 

where Ker (k,M(k);p) is the kernel of the equation. Thus, we need a quadrature rule 
which allows us to approximate 

max 

£AfKer(M4;/>), (32) 
(=l 



where the A,- are the weights of the quadrature and k[ the corresponding nodes, as shown 
in Fig. (10). Here we have used the notation M, = M{k\). In principle, we should be 
interested in the solution of the above expression in the interval p e (k,X). However, 
we can resort to the collocation techniques to solve this equation, imagining that we are 
only interested in the solution M(p) over a finite set of grid points p\,...,pj rather that 
the solution over the complete interval. In such a case, we reach at an equation like 

M{p j )=Y d A i Ker{k i ,M i ;p j ), (33) 
i=\ 

for each pj, with the added advantage that if the set of grid points {pj} is selected as 
the same set of points of the quadrature rule, i.e. pj = kj for j = 1, . . . 7 n max , we are 
left with a system of n max algebraic equations for the n max unknowns M ; . This is the 
feature of collocation techniques, an integral equation over an interval is translated into 
an algebraic system of equations. 

Let us apply these techniques to solve eq. (30) in a convenient interval (?c, A). A good 
integration interval, will be one for which, as K — > and A — > °°, both the infrared and 
ultraviolet behaviors of the mass function we discussed previously have settled over a 
couple of orders of magnitude of momentum. For concreteness, let us select the interval 
(10~ 4 , 10 3 ). The next step is to select a quadrature rule to perform the integration. The 
simplest choice would be a trapezoidal rule and thus a set of grid points in the selected 
interval. A logarithmic grid is the most appropriate, since there are several orders of 
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FIGURE 10. A quadrature rule for integration 



magnitude involved. Furthermore, one has to make sure that the integration points are 
uniformly distributed over the interval of integration. It is desirable to have the same 
number of grid points over each decade in the interval. For example, if we want to 
distribute /i points in each decade of the interval, and we want 

rimax = [logio(A) -log 10 (K)]/i + 1 (34) 

grid points, defining 



we observe that the recursion 

p J+1 = Pj + KeV-W (e 77 - 1 ) , 7 = 1,..., n max - 1 (36) 

provides us with a set of points uniformly distributed over the interval, exactly ji of them 
in each decade. 

We now need to construct the integration weights of our trapezoidal quadrature. 
Notice that each weight A/ in the logarithmic grid should depend of the grid point pj, 



In (10) 
A* 



(35) 



as opposed to the ordinary "equally spaced" weights of a linear grid. Therefore, we can 
define them as 

Ay = ^(pj+i-pj-i) j = 2,...,n max -\ 

1 1 

Ai = ^{P2-Pl), K max = ^{Pn max -p nmax -\) ■ (37) 

With these ingredients, we are ready to write down the system of equations to be solved. 
From the analytical treatment, we expect M(p) to behave like a constant as p — > 0. 
That means that in the denominator of eq. (30), we can write md yn = M\, and thus, for 
j = 1 , . . . , n max , we have an equation of the form 

M J = ^ E Ai T^L \ k Mi- j) + %-*U-i)\, (38) 

since if z > j, then /c; > kj. We can now proceed to solve the system of equations given 
above by standard techniques, like the Newton-Raphson method. 



In order to test the stability of the procedure, we should establish the dependence 
of the solution upon the number of grid points. In Fig. (11) we depict the behavior of 
m^yn as a function of /i, and observe that for /i > 8, the variation of mj yn with the 
number of points per decade is negligible. In Fig. (12), we depict the mass function 
M(p) for different values of /i. We observe precisely that for /i > 10, all the curves 
lay practically on top of each other, which points out the fast stability of the numerical 
procedure. Figure (13) shows de comparison of the analytical expression, eq. (24), and 
the numerical result for ji = 20. The agreement of the curves both in the infrared and 
ultraviolet allows us to readily infer that both the infrared constant behavior of M(p) 
and the ultraviolet fall-off as l/p 2 . The difference, of course, arises in the intermediate 
region, where the effects of the linearization are stronger. 



Although from this exercise we conclude that QED3 is a theory that supports DCSB, 
we still need to find out whether the solution to the SDE under the set of simplifying 
assumptions, eq. (24), supports confinement. Such an issue is discussed below. 



Confinement Test 

Confinement is connected to our impossibility of detecting colored states. From the 
theoretical point of view, there are several criteria connected to the infrared behavior of 
Green's functions to establish confinement (see, for example, [26]), like the Kugo-Ojima 
criterion [27] and the Gribov-Zwanzinger confinement scenario [28], which demand an 
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FIGURE 11. Stability test of the numerical solution with m c i yn 



infrared suppressed behavior of the gluon propagator and quark-gluon vertex and an 
infrared enhancement of the ghost propagator in QCD. Recent lattice simulations, how- 
ever, seem to disagree with those criteria, particularly in the behavior of the ghost and 
gluon propagators [29]. Here, we adopt a different point of view to establish confinement 
in QED3, based upon the Axiom of Reflection Positivity [30], which establishes that if 
an Euclidean Green's function describes a stable excitation, it should be positive defini- 
tive. This statement provides us with a confinement test that would help us understand 
this phenomenon in QED3 and its connection with DCSB [17]. 

First we need to review the physics behind confinement in QED3 [16]. The potential 
between two static charges, as the separation between them, r — > °°, behaves like 
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This means that in the quenched theory, G = 1, the potential is logarithmically con- 
fining. This is precisely the scenario we are considering. Bear in mind, however, that 
the situation changes drastically if loops of massless fermions are taken into account. 
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FIGURE 12. Mass function for various points per decade 



Considering Nf fermion families, we can easily prove that 

STtk 



G(k) 



Snk + e 2 N 







as 







(40) 



which is understandable on physical grounds, since massless fermions cost no energy to 
produce, and thus, they infinitely screen the static charges, sweeping away confinement. 
However, if the fermions within the loops acquire masses either dynamically or explic- 
itly, they provide an effective screening, rendering G(k) ^ and finite as k — > 0, and thus 
confinement is reinstated. 

Confinement can be studied from the fermion propagator in the following form [16]: 
Let us start defining the following space-averaged Schwinger function 



A(t) 



where 
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FIGURE 13. Analytical vs. Numerical comparison of the Mass function 

is the scalar part of the fermion propagator. The Axiom of Reflection Positivity [30] 
establishes that if the propagator describes a stable particle, 

A(f) ~ e-' nt (43) 

and m corresponds to the mass of the particle. Such a behavior is found, for example, 
from the free fermion propagator. Actually, this test is performed in lattice simulations 
to obtain the mass of hadronic resonances. On the other hand, when the propagator 
possesses complex singularities, A(t) is no longer positive definite, but behaves like 

A(f ) ~ e~ mt cos(m 2 t + 8) , (44) 

where m\ and are the position of the said complex poles. In this case, the propagator 
describes a confined excitation. There are alternative formulations of this test [17]. 

If we compute the spatially averaged Schwinger function, eq. (41), with the analytical 
form of our problem, eq. (24), we obtain the curve shown in Fig. (14). We can see that, 
indeed, our solution induces an oscillatory behavior on the Schwinger function, which 
is more evident in Fig. (15), where we plot the logarithm of the absolute value of A(t). 
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FIGURE 14. Confinement test funtion A(f ) 



Each peak corresponds to a solution of A(t) = 0. Thus, we conclude that both DCSB 
and confinement are present in our study. 



A word of caution is at hand. The positiveness of the Green's function is a sufficient 
condition for confinement, but not necessary. The counterexample is QCD in (1+1) di- 
mensions [31]. There, the propagator exhibits complex singularities, but the color sin- 
glet meson bound state amplitudes vanish precisely at the fermion mass shell momenta. 
However, confinement is realized in such a theory via a failure of the so called cluster 
decomposition property [32]. 



FINAL REMARKS 



Through these lectures, we have reviewed the dynamical generation of light quark mass. 
Such a mass is a measure of the confining energy of these quarks inside the nucleons, 
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FIGURE 15. Logarithm of the confinement test funtion A(f ) 



and is determined by the strong interactions of QCD. To understand the relationship 
between confinement and DCSB, we have worked out an neat example of a toy model 
of QCD which exhibits both these phenomena, QED3. The SDE is the platform we have 
selected for this purpose, performing a truncation the tower of SDEs under simplifying 
assumptions. We have presented a detailed analytical solution to the gap equation of 
QED3, and also, we have developed a numerical procedure to look for the dynamically 
generated fermion propagator. A confinement test performed on the propagator revealed 
that this chirally asymmetric solution to the SDE also describes a confined excitation. 
This exercise points out the key features of strong interactions which reveal the origin 
of mass of the visible Universe. 
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